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ABSTRACT 


In  this  report  crosstalk  and  distortion  of  high  speed 
pulses  due  to  multiple  reflections  are  studied.  The  moment 
method  is  used  to  solve  for  the  currents  and  charges  on  a  set 
of  parallel  microstrip  transmission  lines.  Next,  these 
solutions  are  employed  to  construct  the  capacitance  and 
inductance  matrices  for  the  system.  A  modal  analysis  of  an 
n-line  transmission  line  system  is  used  to  evaluate  the  de¬ 
pendence  of  crosstalk  on  line  spacing  and  input  pulse  risetime 
Otner  results  are  given  which  demonstrate  the  difficulty  of 
obtaining  accurate  data  for  poorly  matched  loads. 
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1 .  INTRODUCTION 

This  report  is  an  analysis  of  crosstalk  and  reflections 
in  digital  circuits  on  printed  circuit  (pc)  boards.  Every 
year  packaging  densities  go  up  and  switching  times  come  down 
with  the  result  being  an  increase  in  crosstalk  and  reflections 
in  the  circuits.  Therefore,  a  need  has  risen  for  applying 
microwave  techniques  to  digital  circuits  in  order  to  be  able 
to  predict  crosstalk  and  reflection  voltages  in  a  given 
circuit  which  may  result  in  false  switching  of  the  devices. 

We  begin  in  chapter  two  by  analyzing  a  simple  trans¬ 
mission  line  system  consisting  of  two  parallel  microstrip 
transmission  lines  that  are  infinite  in  length.  The  network 
analysis  of  electromagnetic  fields  is  used  to  derive  a  super¬ 
position  integral  relating  the  charge  distribution  on  a  strip 
to  the  voltage  on  the  strip.  This  is  done  by  finding  a 
Green's  function  due  to  an  arbitrary  current  source  placed 
at  the  z=0  plane  (fig.  1)  (page  5)  in  a  transmission  line  set  up 
in  the  transverse  direction.  Galerkij.'s  method  is  used  to 
find  an  approximate  solution  for  the  charge  distribution. 

The  application  of  Galerkin's  method  resulted  in  several 
integrals  that  could  only  be  done  numerically.  These  integrals 
were  found  to  have  slow  convergence  when  calculated  on  a 
large  Cyber  computer  system.  Asymptotic  approximations 
applied  to  the  integrals  along  with  a  good  choice  of  basis 
functions  resulted  in  fast  convergence.  From  the  calculated 
charge  distribution  the  capacitance  coefficients  are  then 
solved  for. 
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The  same  analysis  is  then  used  to  derive  a  superposition 
integral  relating  the  current  distribution  on  a  strip  to  the 
flux  linking  the  strip  to  the  ground  plane.  Again,  Galerkin's 
method  was  applied  to  give  an  approximate  solution  for  the 
current  distribution  on  the  strip.  The  inductance  coefficients 
are  then  solved  for  from  the  calculated  current  distribution. 

In  chapter  three  a  general  analysis  of  an  N-conductor 
transmission  line  system  is  presented.  The  technique  used 
is  one  of  modal  analysis. 

The  telegrapher  equations  for  a  multi-conductor 
transmission  line  system  are  combined  to  give  an  eigenvalue 
equation.  The  eigenvalues  are  the  inverse  of  the  propagation 
constants  of  the  system  squared  and  are  found  from  the  pro¬ 
duct  of  the  inductance  and  capacitance  matrices  associated 
with  the  system.  The  eigenvectors  represent  the  modal 
voltage  amplitudes.  The  idea  of  modal  transmission  and 
reflection  matrices  is  developed  here. 

The  general  modal  transmission  line  theory  is  then 
specialized  to  the  original  problem  of  two  parallel  micro-* 
strip  transmission  lines  terminated  by  digital  devices.  The 
inductance  and  capacitance  matrices  are  found  through  the 
tecnniques  presented  in  chapter  two. 

A  computer  program  was  written  from  the  results  of 
tne  analysis  given  in  chapters  two  and  three.  The  program 
assumes  tne  input  to  the  system  to  be  a  pulse  with  rise  time 
and  length  specified  by  the  user.  The  program  uses  a  time 


stepping  technique  where  the  input  voltage  ramp  of  the  pulse 
is  approximated  by  a  series  of  step  functions. 

Numerical  results  are  given  for  the  load  voltages  in 
the  system  as  a  function  of  time.  Different  sets  of  results 
are  presented  for  varying  rise  times  of  the  input  pulse  as  well 
as  varying  dimensions  of  the  microstrip  transmission  lines. 


2.  DETERMINATION  OF  CAPACITANCE  AND 


INDUCTANCE  COEFFICIENT  MATRICES 

In  this  chapter,  the  inductance  and  capacitance 
coefficient  matrices  associated  with  a  system  of  two  parallel 
microstrip  transmission  lines  that  are  infinite  in  length 
are  found.  From  elementary  electrostatics  we  have 


[Q] 

=  [Cl  [V] 

(2-la) 

m 

=  [L] [I] 

(2-lb) 

The  problem  can  be  solved  by  the  formulation  of  two 
integral  equations.  One  integral  equation  relates  the  voltage 
on  the  strip  to  the  charge  distribution  on  the  strip,  and 
the  other  relates  the  current  distribution  on  the  strip  to 
the  flux  linking  the  strip  to  the  reference  conductor  ground 
plane.  The  results  from  these  integral  equations  are  then 
used  in  equations  (2-1)  to  solve  for  both  the  capacitance 
and  inductance  matrices. 

2 . 1  Formulation  of  the  Integral  Relating  Voltage  and  Charge. 

Figure  1  shows  the  geometry  of  the  initial  problem. 

The  conductor  is  assumed  to  be  perfectly  conducting  with 
zero  thickness  and  infinite  length.  The  ground  planes  and 
dielectric  are  assumed  to  extend  to  infinity  in  all  directions 
with  the  dielectric  being  lossless.  This  is  a  good  assumption 
for  the  strip  width  being  much  less  than  that  of  the  width 
of  the  ground  planes  and  dielectric. 


The  transverse  fields  in  each  region  are  given  by 


the  following  integral  [1] ,  [8] . 
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Figure  1  -  Microstrip  Cross  Section 


where 


f  2  =  f 


l  x  z0 


K  =  +  y030,  K  =  | k 1 


o,  is  tne  propagation  constant  in  the  y-direction,  xq  ,  yQ,  and  zq 
are  the  x,  y,  and  z  directed  unit  vectors  respectively,  and  2,-1 
and  1*2  represent  TM  waves  (H^-0)  and  TE  waves  ( E 2*0 ) r es pe c t  i ve ly  . 

By  substituting  equation  (2-2)  into  Maxwell's  curl  equations 
we  can  get  the  following  expression 


where 
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Equations  (2-3)  are  the  transmission  line  equations  in  each 
region.  kj ^  is  the  propagation  constant  in  the  z  direction 
for  both  the  TE  and  TM  waves  for  each  media  and  Zp  ^  is 
the  impedance  in  the  z-direction  for  both  the  TE  and  TM 
waves  in  each  media. 

Figure  2  shows  a  transmission  line  in  a  transverse 
section  of  the  microstrip  line.  The  shorts  at  z=d  and 
z— t  correspond  to  the  ground  planes  located  at  those  points. 
The  boundary  conditions  to  be  satisfied  are  as  follows: 


v, (1) 

(d) 

=  0 

( 2-4a) 

(1) 

.  +* 

(2) 

V 

(0  ) 

=  V ,  1 

(0") 

{  2  —  4  b ) 

„  (i) 

.  + 

(2) 

( 2-4c) 
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Boundary  conditions  (2-4a)  and  ( 2  —  4 d )  state  that  the  voltages 


on  the  ground  planes  are  zero.  Boundary  condition  (2-4b) 
states  that  the  voltage  is  continuous  across  the  z=0  plane 
and  boundary  condition  (2-4c)  states  that  the  currents  in 
each  region  are  discontinuous  by  the  test  current,  i,. 

The  voltage  and  current  in  region  1  are  given  by 


V, (1)  (z)  =  Ae 


^  (1) 


-jk .  z  +  ^ 


(1) 


(1) 


(1)  -  1 —  [Ae‘jkf'  2  -  Bejk?'  ’*]. 


I,  ''  (z)  = 


(1) 


(2-5; 

(2-6; 


The  voltage  and  current  in  region  2  are  given  by 


V,(2)(z)  =  Ce*jk^  )z  +  De^k?  2 


(2) 


(z)  =  z  (1)  [Ce 


-jk 


(2) 


2  rv  jk  , 

-  De J 


(2) 


(2-7 

(2-8 


A  negative  sign  in  the  exponent  indicates  a  wave  traveling 
in  the  positive  z  direction. 

By  applying  the  boundary  conditions  (2-4a) ,  ( 2 — 4 b ) , 
(2-4c),  and  (2-4d)  to  equations  (2-5)  through  (2-8)  we  have 
the  following 


Ae 


+  3e 


jk,(I)d 


0, 


A  +  B 


C  +  D, 
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The  test  current  i.  can  be  expressed  by  the  following 
Fourier  integral  in  the  x  direction 

i,  =  -/*f  ,(a,x')  •  I(x')dx'  (2-16) 

i  — ®  t 


where  i (x  ' ) 


x„ix(x)  +  y0iy (x  ') . 


In  digital  circuits  the  bandwidth  of  the  signal  is 
determined  by  the  rise  time  of  the  pulses  being  transmitted. 
For  rise  times  on  the  order  of  one  nsec.,  the  -3  dB  frequency 
is  .35  GHz.  Studies  done  by  Kitazawa  [1]  indicate  that  a 
static  approximation  would  be  quite  useful  at  frequencies 
this  low  on  typical  printed  circuit  boards.  The  static 
approximation  simplifies  the  analysis  a  good  deal. 

Expanding  the  expression  for  the  transverse  H-field 
and  applying  the  approximation  30  -*  0  we  have 


H  (i) 


(i)  ,  _  \  aVo 

=  j  I .  ( a ;  z )  -^-t- 

1  i  a  | 


727 


da 


(2-17) 
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r 

j 

—  CO 


( a ;  z ) 


aX 


da 


The  voltage  on  the  strip  may  be  found  by  either 
integrating  the  E7-  field  in  region  one  from  0  to  d  or  in¬ 
tegrating  the  Ez~field  in  region  two  from  -t  to  0.  The 
choice  here  is  arbitrary  so  the  succeeding  equations  will, 
all  apply  to  region  one  only.  The  field  component  only 
can  be  found  from  the  i=l  TM  field  since  no  Ez  component 
exists  for  the  l=2  TE  field.  Now, 
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S.-oIi  (a;z)  = 


cosh  [  ( z-d)  )  !  a  |  ]  i  ( at) 


T7r  coth  (  |  a  j  t)  +coth  (  |  a  |  d)  }  sinh  (  !  a  |  dT 


i,  (a)  *  f  /Tn  e-^01*  i  (x")dx' 
1  x 


Therefore,  the  transverse  H  field  component  of  the  TM  wave 
is  given  by 


(1)  _  2_1  r®  ,® 


cosh  [  ( 2-d)  lajj 


2t  ^  £rcoth  (  j  a  j  t)  +coth  (  |  a  [ d)  } sinh  (  |  a  |  d) 


ix(x-)eialx'  -  x)dx'da 


The  E 2  component  is  found  from  the  H  field  from 
Maxwell's  curl  equation  which  can  be  expressed  as 


— • -/tt  7  *  x  2„) 

U-  1 


(2-18) 


1  r°°  ,-acosh  [  ( z  -  d)  1  a  1  ]  _ _ _ 

2^  £(1)  '  t  ercoth  (  j  ct  ]  t )  +coth  ( J  a  j  d)  }  sinh  (  j  ct  |  d) 

i  (x-)e-jc,<x-x')dx'do. 


To  get  equation  (2-18)  into  a  form  that  relates  the 
charge  distribution  to  the  electric  field  Ez,  the  integration 
with  respect  to  x"  is  performed  by  parts  and  the  equation 
of  continuity 


j.J(x')  = 


ix(x') 


(2-19) 


.3  used.  Equation  (2-13)  then  becomes 


(1) 


£{/ aT(a;z)e^a(x-X>)dx'i..(x') 


2  TTCU£ 


-  /  ( /aT  (a;  z)  e 


-ia  (x-x ')  . ^ 


dx  ') — ^ — dx'}  da 


where  T  = 


cosh [ ( z-d)  I  a  I ] 


{ ercoth ( | a | t)  +  coth ( | a | d) s inh ( | a j d) 


The  first  term  in  equation  (2-20)  goes  to  zero  since 


(2-20)  J 


i  (  00  )  =  i  (-»)  so  that  equation  (2-20)  becomes 

XX 


Sz 


(i] 


-l 


j  ax 


2-ruue 


(1) 


f  /T-  (a; z) e  ^  -  a  ( -  j  jc  (x  ')  )  dx'da 

-co.  co  1  ja  J 


2tt£ 


jYj —  LoU i  ^a;  a  ^ X  o  (x')  dx'da 


(2-21) 


Now,  the  voltage  on  the  strip  conductor  is  just  the 
integral  of  the  z  component  of  the  electric  field  from  the 
top  ground  plane  to  the  strip. 


V  ( x) 

V  ( x) 


/gEz  (x,  z)  dz  . 


X 


6  (x') 


-ja  (x-x') 


2tt£ 


(1)  ercoth  (  |  a  1 1)  +  coth  (  |  a  |  d)  ;  |  a 


-dadx ' 


(2-22) 


The  structure  of  the  microstrip  used  in  digital  circuits 
such  as  printed  circuit  boards  does  not  have  the  top  ground 
plane  shown  in  figure  1.  It  is  easy  to  transform  the  original 
problem  to  that  of  the  problem  we  are  concerned  with  by  taking 
the  limit  of  V(x)  as  d  goes  to  infinity.  The  only  term  in¬ 
volving  the  constant  d  is  coth(|ujd)  which  goes  to  one  as 
i  goes  to  infinity.  Therefore,  equation  (2-22)  becomes 

•1 


V  f  x) 


X>  .  ,  ^  ^ 

.  G ( a, x , x  )r(x  ) dudx 


(2-23) 


L  3 


I  where 

4 

,  -ja(x-x')  . 

G  ( a,  x|x')  =  — - - - —  (2-24) 

2tte0  e  rcoth  (  |  a  1 1) +1  |  a  | 

1  All  modes  of  propagation  for  the  two  strip  case  can 

be  expressed  as  a  linear  combination  of  the  basis  modes 
called  the  even  mode  and  the  odd  mode.  The  even  mode  is  the 
mode  for  which  the  potential  on  the  two  strips  is  equal. 

The  odd  mode  is  the  mode  for  which  the  potential  on  the  two 
strips  is  equal  in  amplitude  but  opnosite  in  sign.  There¬ 
fore,  symmetry  can  be  used  to  simplify  the  expression  for 
voltage.  Equations  (2-23)  and  (2-24)  then  become 


V(x)  =  /™/G  (a,x|x'’)a(x")  dadx ' 


(2-25) 


where 


G  (a,  x  |  x' 


2  F  (a) 


cos  ix  cosax'  (2-26) 

for  the  even  mode 


2  F  (a) 


„  ct 


sinax  sinax 


for  the  odd  mode 


and 


trcoth(at)+l 


Since  charge  will  only  exist  on  the  conductor,  then 
the  integration  with  respect  to  x "  is  zero  at  all  points  on 
the  plane  not  occupied  by  the  conductor.  Equation  (2-25)  then 


becomes 


(2-27) 


V(x)  =  /°/0G (a,x | x") a (x ' j dadx" 
where  G  is  as  given  before. 

2 . 2  Solution  of  the  Voltage  Charge  Integral  Equation. 

Prior  to  the  advent  of  large  computers,  engineers  and 
scientists  would  try  to  find  some  form  of  a  given  integral 
equation  which  could  be  solved  analytically.  Now,  computers 
can  be  used  to  find  highly  accurate  "approximate"  solutions. 
The  integral  equation  that  was  derived  in  the  previous  section 
relating  voltage  to  charge  would  be  very  difficult  if  not 
impossible  to  solve.  The  method  that  is  commonly  used  to 
solve  these  problems  is  called  the  method  of  moments. 

It  is  assumed  that  the  reader  has  some  background 
in  linear  spaces.  A  more  detailed  explanation  is  given  in 
reference  [7],  Given, 

L(u)  *  s 

where  s  is  known  and  u  is  unknown,  then  u  may  be  found  by 
a  =  L"1(s) . 

To  do  this,  i-  is  expanded  in  a  set  of  linearly  in¬ 
dependent  basis  functions 
N 


Since  any  function  may  be  exactly  represented  as  an  infinite 
series,  then  the  finite  expansion  for  u  is  only  an  approximation 


A  set  of  testing  functions  is  selected  in  the  domain 
of  L,  The  inner  product  of  L{  u)  *  s  is  evaluated  for  each  wm 
Therefore/  we  have 
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n=l  n 
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The  resulting  set  of  linear  algebraic  equations  can  then  be 
solved  for  the  unknown  coefficients. 

The  method  used  here  is  a  special  case  of  the  method 
of  moments.  The  testing  functions  selected  are  the  same 
as  the  basis  functions.  In  this  case  the  method  of  mements 
is  properly  called  Galerkin's  method. 

The  unknown  function  in  the  case  of  equation  (2-27)  i 
the  charge  distribution  a(x).  The  exact  solution  is  given 
for  a  microstrip  in  free  space  over  a  ground  plane  as 


A  is  a  constant  determined  by  the  potential  of  the  strip  with 
respect  to  the  ground  plane.  For  the  perturbed  strip  over  a 
dielectric  substrate  in  the  presence  of  another  strip,  vari¬ 
ational  terms  can  be  added  to  give  a  good  approximation  of 
the  charge  distribution. 


The  basis  functions  selected  in  this  case  are 


2  (x-s)  ,  2 


'2  -"F(a)  ,b(x*-s) slnax* 


dx 


-b  (x-s ) 
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'  ^  , 2  ? X-3)  | 2 
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Several  of  these  integrals  can  be  analytically  integrated 
and  are  given  as  follows 
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(2-34) 


x-s 
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(2-35) 


rb  cos  ax^ _ 


cos(as) J0(2p)  =  P2  (2-36) 


^.b(x-s)  cosax  dx 

a  A"  r2(x-s)  72 
1  w  J 


sin(as)  J.  (^)  =  P,  (2-37) 
4  12  3 


rb  sin  ax 


ttw  .  ,  ,  T  -aw, 

■y  sin(as)J0  (y) 


(2-38) 


4b(x-s)sinax  dx 

/.  -2 (x-s) ,2 

1  w 


ttw  ,  ,  T  -  wa , 

-y cos  (as)  J1  (y) 


(2-39) 


From  equations  (2-30)  through  (2-33)  and  equations 
(2-34)  through  (2-39)  we  have 


2J±  rtSSLfii*  +  ?!i  /*£I24  .  do 


p -  =  ~ j0  :-p, 

1  tt£  „  °  a  2 


- JT  -)  XT  ~  V. 

a  3  2 


(2-40) 


0  ,  m.,  dc,  ♦  !!l  r»£i*4  2da 
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-P  =  P  2da  +  ^  /“f (a) 

P1  Te.°  a  P4  aa  ttSo  :^Tlp5P4dci 

o  =  .fEiilpp,  da  +  ^  ;*Ei“Lp  2da 

ttg  0  0  a  45  *rc  0  0  a  5 
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(2-43) 


Six  different  integrals  must  be  solved  in  this  set 
of  equations  in  order  to  solve  for  the  unknown  coefficients. 

The  integrals  were  first  evaluated  directly  using  the  IMSL  com¬ 
puter  subroutine  DCADRE.  The  cost  of  evaluating  these  integrals 


was  too  high  because  of  slow  convergence  for  each  integral. 
Asymptotic  approximations  of  the  six  functions  to  be  in¬ 
tegrated  can  be  found  in  a  table  of  integrals  [9]  and  are 
given  as  follows: 

J.-F1C0  p^2^a  , 

ttw  , °°  cos2as  sin  wa  +  cos  2as  +  sin(wa)  +  1  . 

- TT7 - 

-  HiTTIT  e'23e  +  HIe'WSgi("6)  - 

+  51  +  e"wS{e2s3Ei[6(w-2s) ]  +  e“2 S8Ei [ 3 (w+2s) ] } 

-  ^  ew6{e2s8Ei[-S(w+2s) ]  +  e_2s8Ei [ S ( 2s-w) ] }  }  (2-44) 
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In  each  equation,  3  is  a  constant  that  was  added  to  the 
denominator  so  that  each  integral  could  be  calculated  analytically 
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The  accuracy  of  each  equation  is  therefore  dependent  on 
small  8.  The  method  to  calculate  each  integral  is  given  as 
follows . 

Let  A (a)  be  the  asymptotic  approximation  of  the 
function  f{a)  to  be  integrated.  Then  the  integral  can  be 
approximated  as 

I  =  J^af(ct)-  A(a)da  +  A(a)da. 


Ie  =  -4-  L  ^4^P,Pcda 
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Then  we  have  for  the  even  mode 


TTW 

2 

*1 

M 

^  1 

A1 

0 

— 

}2 

I3_ 

— -j, 

_ l 

and  for  the  odd  mode 
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These  equations  are  easily  solved  by  either  matrix  inversion 
or  algebraic  manipulation. 

2.3  Caoacitance  Coefficients 


The  total  charge  on  the  conductor  is  given  as 
Q  =  /b  a(x)dx. 

d 


This  integral  is  easily  evaluated  for  both  the  even  and  odd 
modes.  We  have  for  the  even  mode 
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and  similarly  for  the  odd  mode 


Q'  '  =  A^-yO  +  B2{0) 


(2-53) 


In  a  multiconductor  transmission  line  with  transverse 


fields  [10J 


Q.  =  Z  C. .V. 

1  j-1  «  3 


where  are  Maxwell's  capacitance  coefficients.  In  our 
case  n=2  so  that  we  have  for  the  even  mode 


C11  C12  1 


C  C  1 

21  22  1 
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and  for  the  odd  mode 


Cll-  C12 


C21  C22 


(2-55) 


Using  the  symmetry  of  the  problem  the  charges  on  the  two 
strips  are  related  by 


0,'1’  .  Q2a> 


Q  „  _q  121 


Equations  (2-54)  and  (2-55)  together  with  these  cond 


i tions 


can  then  be  used  to  solve  for  the  capacitance  coefficients. 
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2.4  Formulation  of  the  Integral  Relating  Inductance  and  Current 


The  formulation  here  is  much  the  same  as  that  given 
for  the  charge  voltage  integral  equation.  Figure  1  shows 
the  geometry  of  the  problem.  From  equation  (2-2)  the  trans¬ 
verse  H  fields  are  given  as 


„( 2)  T  (2)  ra(-flB)  +  y(a)  ,  j 

Ht  J-oo  I1  f  -  .  . 

J a2  I  62  m 


.-Jcuw3  s*y 


+  /”  I<21  j.  e-j“xdae'jB"y 

2  ✓  2  72“  2TT 

a  +  S 

where  l|2^  and  I^^  are  as  given  in  equation  (2-15)  . 

Since  it  is  desired  to  find  the  flux  linking  the  ground 
plane  to  the  strip,  then  only  the  x  component  of  the  H  field 
is  needed.  Applying  the  zero  frequency  approximation  we  have 


Kx'21  -  £=■ 

z  /2tt  k  ■ 


(2-56! 


where 


cosh  (  |  a  i  ( t+  z)  )  i  (a) 
cosh  (  j  i ;  t)  +  sinh  (  a  j  t)  coth  (  HiiTd] 
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The  flux  linking  the  strip  to  the  ground  plane  is 
<t>  (x)  =  /°t  U0Hxdz 
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Using  the  symmetry  of  the  problem  we  finally  have 
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where 
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for  the  even  mode  and 

„  ,  . ,  2-0  F  (  a ) 

G(a,x  x  )  =  — “■  -  sin  ax  sin  ax 

t  a 


for  the  odd  mode  and 


That  tr.  j  equation  is  very  similar  to  the  one  given 
for  voltage  an  tr.e  strip  should  not  be  surprising  since  the 
two  problems  art  : jj 1  to  each  other.  Therefore,  the  same 
technique  use!  to  solve  for  the  charge  distribution  on  one 
strip 


due  to  a  voltage  on  the  strip  can  also  be  used  to  solve 
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for  the  current  distribution  on  the  strip  due  to  a  flux 
linking  the  strip  with  ground. 

2 . 5  Inductance  Coefficients 

To  find  the  current  distribution  on  the  strip,  the 
same  method  is  used  as  that  used  for  finding  the  charge 
distribution.  Once  the  current  distribution  is  known  then 
the  total  current  is  found  by 

I  =  /b  i  (x)dx. 
a  y 

In  an  N  conductor  transmission  system,  the  magnetic 
flux  linking  the  n—  conductor  and  ground  is  the  resultant 
of  currents  flowing  in  the  N  independent  conductors  [10] 


where  the  L  are  the  line  coefficients  of  inductance .  In 
nm 

finding  the  currents  on  a  strip,  the  magnetic  flux  was  assumed 
to  have  a  certain  value 


for  the  even  case  and 


5l  =  1,  >2  =  -1 

for  the  odd  case.  Therefore,  for  the  even  case 
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and  for  the  odd  case 


1 

L11 

L12 

-1 

_L21 

L22_ 

^2). 

Again  using  the  symmetry  of  the  problem,  the  currents  on 
the  strips  are  related  by 


I 


1 


(1) 


(1) 


for  the  even  case  and 


for  the  odd  case.  Equations  (2-60)  and  (2-61)  together  with 
these  conditions  can  then  be  solved  for  the  inductance  co- 
ef  f  icients . 


3.  D.  C.  TRANSIENT  ANALYSIS  OF  COUPLED  LOSSLESS 
TRANSMISSION  LINES 

The  D.  C.  transient  analysis  of  coupled  microstrip 
transmission  lines  is  covered  in  detail  in  the  literature.  In 
this  thesis  an  effort  has  been  made  to  apply  these  techniques 
to  real  world  digital  circuits.  The  loads  are  the  digital 
devices  and  the  transmission  line  is  the  microstrip  conductor 
connecting  two  devices  on  a  printed  circuit  board.  The  output 
and  input  impedances  of  digital  circuits  are  highly  non-linear 
and  the  reflections  from  these  types  of  loads  can  cause  dis¬ 
tortion  to  the  original  signal  so  that  false  switching  of  the 
receiving-device  may  occur.  In  addition,  coupling  to  adjacent  lines 
occurs  from  the  drive-line  which  can  cause  false  switching  to 
occur  on  the  non-drive-lines.  In  this  chapter  a  method  is  pre¬ 
sented  which  can  be  used  to  calculate  distortion  from  both 
coupling  and  reflections. 

3 . 1  Eigenmode  Expansion  Method  for  Multi-Conductor 
Transmission  Lines 

Several  papers  [2,3,4]  were  published  during  the  60' s 
and  early  70' s  on  crosstalk  between  conductors  in  three  wire 
transmission  line  systems.  The  models  used  are  not  easily 
extendable  to  more  than  three  conductor  transmission  lines  and 
do  not  take  into  account  reflections  from  the  terminating 
impedances.  Recently,  some  techniques  using  eigenmode  expansion 
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methods  have  been  presented  [5,6].  The  advantage  of  a  modal 
analysis  technique  is  that  an  N  conductor  transmission  line 
system  can  be  modeled  as  a  single  tube  terminated  by  a  network 
through  which  N  modes  of  propagation  are  possible. 

We  begin  by  defining  the  following  matrices  which  give 
the  line  voltages  and  currents  of  the  system  (fig.  3)  and  the 
capacitance  and  inductance  coefficients  coupling  those  lines. 
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(3-1) 


(3-2) 


(3-3) 


(3-4) 


The  telegrapher  equations  for  an  N+l  conductor  trans 
mission  line  system  are  given  as 


^(V(y,t)]  =  -[Lj^[I(y,t)]  ( 3-5a) 

! 

^(I(y,t)]  =  -[C]^-[V(y,t)  ]  . 

l 


( 3-5b) 


ciurc  3  -  Multiconductor  Transmission  Line  System 
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Equations  (3-5)  can  then  be  combined  to  give  the  wave  equations 
for  both  voltage  and  current. 


2  2 
-^[V(y,t)]  =  [L]  tC]-2-=-[V(y,t)] 
ay  at 


(3-6a) 


2  2 
■^IKy/t)]  =  1C]  [L]-2-r[I(y,t)l 
ay*  at 


(3-6b) 


It  should  be  stressed  that  these  equations  are  for  a  TEM  mode  of 
propagation  only.  At  high  frequencies  the  mode  of  propagation 
is  not  TEM  for  a  microstrip  structure.  However,  at  low  fre¬ 
quencies  the  propagation  is  very  nearly  TEM  so  that  very  good 
calculations  can  be  gained  with  this  method. 

Now,  define  matrices  for  generator  voltages  and  source 
and  load  imnedances. 
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(3-8) 


[vs (t) ]  = 


'SI 
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SN 


(3-9) 


The  source  and  load  matrices  describe  the  networks  that  termi¬ 
nate  each  end  of  the  system  (fig.  3)  ,  and  the  source  matrix 
gives  a  description  of  all  the  generators  at  the  y=0  end  of  the 
system. 

Solutions  to  (3-6)  are  given  as 

fV (y / 1)  ]  -  [VQ]  •  g(t  ±  £>  (3-10a) 

[I (y»t) ]  =  [IQ]  •  g(t  ±  £)  (3-10b) 

The  variable  v  is  the  wave  velocity.  A  positive  sign  in  front 
of  ^  denotes  a  wave  propagating  in  the  negative  y  direction  and 
a  negative  sign  denotes  a  wave  propagating  in  the  positive  y 
direction.  low  we  substitute  (3-10a)  into  (3-6a)  to  get 

-~[VQ]  g  (t  t  £)  =  [L] [C] [VQ]  g  (t  ±  £).  (3-11) 

Equation  (3-11)  is  an  eigenvalue  equation  where  — are  the  eigen 
.  .  v 

values  and  [VQ  ]  are  the  eigenvectors.  Since  for  an  N  con¬ 
ductor  transmission  line  system  [L] [C]  is  an  NxN  matrix,  then 
there  will  be  N  eigenvalues  (not  necessarily  distinct)  and  N 
eigenvectors.  The  eigenvalues  and  eigenvectors  are  related  to 

the  modes  of  the  system,  each  mode  having  a  velocity  —  and 

(in )  ^ 

amplitude  [V  '  ] .  Thus  we  can  model  our  system  of  N  conductors 

as  a  single  tube  with  N  modes  of  propagation.  Similarly,  from 

equations  (3-10b)  and  (3-6b)  an  eigenvalue  equation  for  the 
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currents  can  be  found  with  — ~  again  being  the  eigenvalues  and 
(m)  v  m 

[I  ]  being  the  eigenvectors  corresponding  to  the  eigenvalues. 
Let 

[Sv]  =  {  CV^1)  ]  ,  [V^2)],  [V<N>]}  (3-12) 

and 

[S_]  =  {[I*1*],  [I^2)  ]  ,  [I,!N)]}  •  (3-13) 

loo  o 

From  equation  (3-5a)  the  modal  voltages  and  currents  are  related 
by 

tO  =  ~  CL]  [Iom] 
m 

or  from  (3-12)  and  (3-13) 

[SI]  =  [L]_1(SV] [A]  (3-14) 

where  [A]  is  a  diagonal  matrix  whose  elements  are  the  eigen¬ 
values  — —  of  the  matrix  [L]  [C]  . 
vm 

The  total  voltage  at  any  point  on  the  line  is  the  Siam  of 
the  incident  and  reflected  voltages 

CV(y,t)J  =  Cvinc(y,t)]  +  (vref(y,t)]  . 

Thus,  the  modal  voltages  can  be  related  to  the  conductor  volt¬ 
ages  by 
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Since  the  generator  voltages  are  located  at  y=o,  then  we  can 
write 


tVinc(o't)3  =  [ts]  tvs(o't)]  +  tps]  tVref (o't}]  (3"19) 

where  t  and  p  are  the  transmission  and  reflection  coefficients 
s  s 

at  the  source  end,  respectively.  [x  ]  and  [p  ]  are  defined  by 

S  b 

[Tg]  =  [ZClURg]  +  [ZC]}"1  (3-20) 

and 

[pj  =  {[R  ]  -  ( ZC] } • {  [R  ]  +  [ZC]}"1.  (3-21) 

S  S  2 

From  (3-18)  and  (3-19) 

[Sv][Finc(o't)3  *  CTsHVs(o,t)]  +  tps]tVref(o't)] 
or 

[Finc(o't)]  =  tsvr1(TS]  [Vs(o't)]  +  [Sv]’1[ps]  * 

I*ref  (o't>) 

=  [Sv]"1[is] [Vg(ort)]  +  [Sv]_1[ps]  • 

[Sv][Fref].  (3-22) 

Let 


(3-24) 


IPs . 3  =  fSv3  fps3[Sv3 


be  the  modal  reflection  and  transmission  coefficients  at  the 
source  end.  Then  (3-22)  can  be  written 


[Finc(°'t)1  *  fTs(m)]IVa(o,t))  +  [Ps<m)][Fref(o,t)J 


(3-25) 


Similarly,  at  the  load  end  of  the  line 

[V(l,t)]  =  [Sv]([Finc(l,t)]  +  [Fref  (l,t)] }.  (3-26) 

Now 


tVref(1,t)1  =  [pL] [Vinc(1't)]  (3-27) 

where  [p^]  is  the  voltage  reflection  coefficient  matrix  at  the 
load  end  of  the  line  given  by 


[PL1  =  {  [Rl1  -  [ZC]  }  •  (  [FL^]  +  [ZC]  } 


-1 


(3-28) 


From  (3-26)  and  (3-27) 


lsv11Fref(1't)1  -  [0L1[vinc(1't)1 


or 


[F _ e  (1  ,  t)  ]  =  [S.J"1(PT1  [V,__a  r  t)  1  • 


(3-29) 
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Thus,  we  can  define  the  modal  reflection  coefficient  matrix  at 
the  load  end  of  the  system  by 

[pL(m)]  =  [Sv] "1 f  pL] .  (3-30) 

Therefore, 

tFref(1't,]  *  [PL(m)]  fFinc(1't)1  *  (3-31) 

It  can  be  seen  that  once  the  source  voltage  matrix,  the 
terminating  impedances  and  the  inductance  and  capacitance 
coefficient  matrices  are  known,  then  we  can  find  the  voltages  on 
the  lines  in  our  system. 


1 


4 
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4 .  RESULTS 

From  the  equations  developed  in  chapters  2  and  3,  a 
program  was  written  to  evaluate  the  voltages  on  both  ends 
of  a  system  of  two  Darallel  microstrio  transmission  lines 
with  any  type  of  resistive  loads.  This  chanter  gives  some 
results  from  that  program  as  well  as  some  experimental 
results . 

4.1  Assumptions  and  Simplifications. 

We  would  like  to  be  able  to  find  the  response  of  a 
system  to  a  pulse  with  finite  fall  and  rise  times  (figure 
4).  For  computational  purposes  an  approximation  to  this 
pulse  can  be  used  where  the  voltage  ramps  at  the  beainnina 
and  end  of  the  pulse  can  be  approximated  by  a  series  of  Dosi 
tive  and  negative  going  voltacre  steos,  each  with  an  amplitud 
and  length  determined  by  the  rise  time  and  amplitude  of  the 
original  pulse  (figure  5).  The  time  between  steos,  it,  is 
given  by  the  rise  time  divided  by  some  constant  c.  Simi¬ 
larly,  the  amplitude,  Av,  of  each  step  is  given  by  the  ampli 
tude»  VQ>  of  the  original  pulse  divided  by  the  same  constant 
The  larger  the  constant  c,  the  better  the  approximation 
becomes . 

Therefore,  our  problem  essentially  becomes  one  of 
exciting  the  transmission  line  system  with  a  series  of  step 


functions.  From  equation  3-25,  it  can  be  seen  that  both 
modes  of  the  system  will  be  excited  by  a  unit  step  gener¬ 
ated  on  only  one  line.  These  model  voltages  can  be  viewed 
as  voltage  discontinuities  proDagating  down  a  tube  with 
different  velocities. 

Consider  mode  1  of  our  system  only.  At  time  t  = 
1/v^,  mode  1  will  be  incident  upon  the  far  end  load.  Ex¬ 
panding  3-31  for  the  two  line  system. 


—  “ 

^  “ 

F2ref  (1 '  1/vl) 

(m) 

°L11 

A  ^ 

PL12 

F1inc  (1  ,l/v1 ) 

F2ref  (1 '  1/vi) 

(m) 

°L21 

A  ^ 

PL22 

0 

it  can  be  seen  that  both  modes  will  be  reflected  back  even 
though  only  one  mode  has  been  incident  on  the  far  end  load. 
After  both  modal  voltage  discontinuities  have  been  incident 
on  the  far  end  load,  four  voltage  discontinuities  will  be 
prooaqating  back  towards  the  source  end  of  the  tube.  Two 
voltage  discontinuities  will  be  in  mode  one  and  the  other 
two  in  mode  two . 

In  general,  for  an  N  conductor  transmission  line  sys¬ 
tem  excited  by  a  unit  step  on  one  line,  the  number  of  voltage 
discontinuities  all  prooagatinq  in  the  same  direction  down 
a  tube  for  a  single  mode  is 
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where 

N  *  number  of  voltage  discontinuities  existing 
m  along  a  tube  for  the  mth  mode 

N  ^  =  number  of  reflections  that  have  occurred 
*  for  all  the  modes. 

It  can  be  seen  that  as  time  goes  on,  an  expanding 
amount  of  computer  memory  will  be  needed  for  the  bookkeeDing- 
so  to  speak.  In  the  current  version  of  the  program  none  of 
the  cross-coupled  modes  are  propagated  back  down  the  tube 
past  the  second  reflection.  The  results  to  be  nresented  show 
where  this  aoproximation  may  or  mav  not  work. 

4 . 2  Brief  Program  Description 

The  program  consists  of  four  basic  parts.  The  first 
oart  of  the  proaram  provides  interaction  with  a  user.  The 
user  is  promoted  for  data  concerning  the  geometry  of  the 
structure,  the  input  pulse  width  and  rise  time,  and  data 
describing  the  loads  at  the  ends  of  the  lines. 

The  second  part  of  the  proaram  calculates  the  induc¬ 
tance  and  capacitance  coefficient  matrices  using  the  eauations 
derived  in  chapter  two.  This  part  of  the  proaram  takes  more 
cou  time  than  the  rest  of  the  program  combined  due  to  the 
integrals  (eauations  2-40,  2-41,  2-42,  and  2-43).  This  • 
oart  of  the  program  is  not  easily  adaptable  to  microstrip 
transmission  line  systems  with  more  than  two  signal  conductors. 


This  is  due  to  the  increasing  complexity  and  number  of  inte¬ 
grals  that  would  have  to  be  computed. 

The  third  part  of  the  procrram  computes  the  eigen¬ 
values  and  eigenvectors  of  the  product  of  the  inductance 
and  capacitance  coefficient  matrices.  Initial  values  for  the 
modal  transmission  and  reflection  matrices  are  also  computed 
here  assuming  that  the  system  has  not  yet  been  excited. 

The  fourth  part  of  the  program  consists  of  a  loop 
which  generates  voltage  steps  into  the  system  and  propaqates 
the  voltage  steps  through  the  system.  Each  iteration  through 
the  loop  represents  an  increment  of  time  At  given  in  figure  5. 
At  the  end  of  each  loop  new  modal  transmission  and  reflection 
matrices  are  calculated  and  the  line  load  voltaaes  are  calcu¬ 
lated  from  the  modal  load  voltaaes  and  stored.  None  of  the 
cross-couoled  modes  after  the  second  reflection  are  propa¬ 
gated. 

The  third  and  fourth  parts  of  the  program  are  easily 
adaptable  to  transmission  line  systems  with  more  than  two 
signal  conductors. 

4 . 3  Linear  Loads 

Consider  figure  6.  The  transmission  lines  are  micro¬ 
strip  lines.  The  results  in  fiqures  7-9  all  have  the  follow¬ 


ing  parameters: 


R2  -  R3  -  R4  *  100  0 


er  =  5.0 


1  «  .3  m 

h  *  w  *  10  mils 
pulse  width  =  6.0  ns. 


The  results  in  figures  7  and  8  show  the  effects  of 
increasing  the  separation  between  strips.  The  input  pulse 
has  rise  and  fall  times  equal  to  one  nanosecond.  Note  the 
markedly  decreased  induced  voltages  on  the  passive  line  when 
the  spacing  is  doubled.  The  crosstalk  has  fallen  by  about  a 
factor  of  two. 

The  results  in  figures  7  and  9  show  the  effects  of 
shortening  the  rise  and  fall  times  of  the  inout  Dulse.  The 
spacing  in  both  figures  7  and  9  is  10  mils.  The  rise  and 
fall  times  of  the  input  pulse  in  figure  7  is  one  nanosecond 
while  the  rise  and  fall  times  of  the  inout  pulse  in  figure  9 
is  .5  nanoseconds.  The  crosstalk  increases  substantially 
with  decreasing  rise  and  fall  times  showing  the  dependence 
of  the  crosstalk  on  the  derivative  of  the  input  pulse. 

In  all  cases  the  first-order  reflections  occurri no 
at  about  3.5  nanoseconds  and  about  10  nanoseconds  have  caused 
substantial  shifts  in  the  load  voltaaes.  Any  later  reflec¬ 
tions  are  well  damped.  Thus,  we  can  say  that  in  these  cases 
the  simplification  of  not  propagating  any  of  the 
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Figure  7  -  Time  Response  of  Two  Line  Transmission  Line  System 
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Figure  8  -  Time  Response  of  Two  Line  Transmission  Line  System  With  Increased  Spacing 
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Response  of  Two  Line  Transmission  Line  System  With  Decreased  Rise  Time 


cross-coupled  modes  after  the  second  reflection  will  not 
significantly  alter  the  results. 

4 . 4  Experimental  and  Theoretical  Results 

For  figures  10  and  11  the  following  parameters  were 

used : 

R1  •  50  0 

R2  =  R3  =  R4  =  88  Q 
er  =  5.0 
1  =  .381  m 

h  =21  mils 
w  =20  mils 
S  =5  mils 

rise  time  =  1.0  nanoseconds 
input  voltage  =  .4  volts 

The  pulse  width  for  the  theoretical  data  shown  in 
figure  10  is  only  15  ns  as  compared  with  that  for  the  experi¬ 
mental  results  which  used  an  input  pulse  width  of  about  30  ns 
This  was  done  because  of  the  limitation  of  the  number  of 
points  that  could  be  plotted  by  the  plot  routine  being  used. 

Other  than  the  difference  in  pulse  widths,  the  theo¬ 
retical  data  correspond  very  well  with  the  experimental  . 
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Figure  10  -  Theoretical  Results  for  Transmission  Line  System  With  Well  Matched  Loads 
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Figure  11  -  Experimental  Results  for  Transmission  Line 
System  With  Well  Matched  Loads 


Figure  12  -  Experimental  Results  for  Transmission  Line 
System  With  Poorly  Matched  Loads 


Figure  13  -  Theoretical  Results  for  Transmission  Line  System  With  Poorly  Matched  Loads 
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This  indicates  that  the  input  pulse  was  well  damped 
due  to  loads  that  were  fairly  well  matched  to  the  line 
impedances . 

In  figures  12  and  13  the  following  parameters  were 

used: 

R1  =  50  a 

R2  =  R3  =  R4  =  open  circuit 
cr  =  5.0 
1  =  .381  m 

h  =21  mils 
w  =20  mils 
s  =10  mils 

rise  time  =  1.0  nanoseconds 
pulse  width  =  30  nanoseconds 
input  voltage  =  .4  volts. 

Here  we  see  a  breakdown  in  the  theoretical  results 
past  about  10  ns.  This  was  due  to  not  propagating  the  cross 
coupled  modal  voltages  after  the  second  reflection.  The 
loads  are  poorly  matched  and  severe  ringing  occurs  in  this 


case. 


5.  CONCLUSIONS 


In  this  report,  we  nave  studied  the  problem  of  cal¬ 
culating  voltages  at  the  ends  of  a  three-conductor  micro¬ 
strip  transmission  line  system.  The  problem  of  calculating 
the  capacitance  and  inductance  coefficient  matrices  was  first 
solved.  We  first  derived  an  integral  equation  relating  vol¬ 
tage  and  charge  by  analyzing  a  simple  transmission  line  model 
3et  up  in  the  transverse  direction  of  a  microstrip  structure. 

In  a  similar  fashion,  an  integral  eauation  relatincr  magnetic 
flux  and  current  was  derived.  The  method  of  moments  was 
used  to  solve  for  the  unknown  charges  and  currents.  The 
resulting  integrals  were  found  to  be  very  expensive  to  inte¬ 
grate  so  asymptotic  methods  were  used  to  shorten  the  amount 
of  computer  time  needed  to  solve  these  integrals.  Once  the 
currents  and  charges  were  solved  for,  then  the  capacitance 
and  inductance  coefficient  matrices  were  easily  solved  for. 

A  method  for  findina  the  voltages  at  the  loads  of 
an  N  conductor  transmission  line  system  was  presented.  The 
telegrapher  equations  were  combined  to  yield  a  second-order 
partial  differential  equation.  A  general  solution  to  this 
equation  was  then  substituted  into  the  differential  equation 
yielding  an  eigenvalue  equation  with  the  principal  matrix 
being  the  product  of  the  inductance  and  capacitance  coefficient 


matrices  associated  with  the  system.  The  resulting  eigen¬ 
values  were  one  over  the  propagation  velocities  squared, 
and  the  eigenvectors  were  the  modal  amplitudes  of  the  system 
The  result  was  that  an  N-conductor,  transmission  line  system 
could  be  modeled  as  a  single  tube  with  N  modes  of  propaga¬ 
tion.  The  idea  of  modal  transmission  and  reflection  coeffi¬ 
cient  matrices  was  developed.  The  modal  load  vol tapes  were 
then  solved  for  by  relating  the  incident  and  reflected  modal 
voltages  throuah  the  modal  transmission  and  coefficient 
matrices . 

In  chapter  four  some  results  were  presented  which 
showed  some  of  the  difficulties  of  using  time-domain  methods 
for  analyzing  a  transmission  line  system.  The  chief  problem 
was  the  need  for  an  expanding  memory  to  keep  track  of  cross - 
coupleu  inodes  produced  when  a  single  mode  was  incident  upon 
a  load.  For  loads  fairly  well  matched  to  a  transmission 
line  system,  an  input  pulse  is  damped  quickly  enough  so  that 
cross-coupled  modes  need  not  be  propagated  after  the  second 
reflection . 
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